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We consider the modified Ising model introduced by de Oliveira et al. [J.Phys.A 26, 2317 (1993)], 
where the temperature depends locally on the spin configuration and detailed balance and local 
equilibrium are not obeyed. We derive a relation between the linear response function and correlation 
functions which generalizes the fluctuation-dissipation theorem. In the stationary states of the 
model, which are the counterparts of the Ising equilibrium states, the fluctuation-dissipation theorem 
breaks down due to the lack of time reversal invariance. In the non-stationary phase ordering kinetics 
the parametric plot of the integrated response function x(t>tw) versus the autocorrelation function 
is different from that of the kinetic Ising model. However, splitting x(t>tw) into a stationary and 
an aging term x(t, tw) = Xst(t — t w ) + Xag(t, t w ), we find Xag(t, ~ t w ax f(t/t w ), and a numerical 
value of a x consistent with a x = 1/4, as in the kinetic Ising model. 

PACS: 05.70.Ln, 75.40.Gb, 05.40.-a 

I. INTRODUCTION 

In equilibrium systems the integrated response function \(t,t w ) and the autocorrelation function C(t, tm) depend 
only on the two time difference r = t — t w and are related by the fluctuation dissipation theorem (FDT) [l| 

x(r)=2(C), (1) 

where 

x(C)=P[C(0)-C(r)] (2) 

and (3 = —dx(C)/dC is the inverse equilibrium temperature. 

In recent times many studies have considered the possibility to relate response and correlation functions in non 
equilibrium systems. In aging systems, such as glassy materials and coarsening systems, relaxation properties depend 
both on t w and t and FDT breaks down. In this context, guided by the solution of mean field spin glass models, 
Cugliandolo and Kurchan proposed that, despite the explicit two-time dependence of X an d C, a relation analogous 
to Eq. QJ, namely 

X (t,t w )=x(C), (3) 

may still hold for a large class of systems in the large t w limit. The functional form of x(C), however, is different 
from the equilibrium one @ and system dependent. In particular, it was shown 3l that x(C) is related to basic 
properties of the equilibrium states; because of that, aging systems can be classified Q into few classes according to 
O ■ the shape of x(C). Phase ordering systems are characterized by a broken line shape of x(C). More precisely, for C 
larger than the Edwards- Anderson order parameter qEA, FDT still holds. For C < qEA one has an horizontal line, 
^ i namely a constant integrated response function. In analogy to equilibrium systems, the quantity —dx{C)/dC can be 
interpreted Q as an effective inverse temperature /3 e //(C) of the non equilibrium state. In coarsening systems this 
quantity takes two values, the temperature of the reservoir /3 e //(C) = 0, in the region C > qEA, and /3 e //(C) = 0, 
for C < qEA- The feature /3 e //(C) — (3 in the region of the largest values of C(t, t w ) is quite general in systems where 
local equilibrium |6( is obeyed. Local equilibrium, in fact, implies that on short timescales FDT holds. Since small 
time separations corresponds to the largest values of C(t,t w ), in this regime one has (3 e ff(C) = (3. The situation is 
different in systems where local equilibrium is not obeyed. In this case one does not expect to observe FDT even in 
the short timescale regime and the definition of a thermodynamic temperature from %(C) may be incorrect. 

In this paper we study the fluctuation dissipation relation in a two dimensional Ising model without detailed balance 
(IWDB), originally introduced in Ref. 0, for which local equilibrium does not hold. This spin system is analogous to 
the kinetic Ising model (KIM) but the temperature entering the transition rates depends on space and time through 
the system configuration. This model is known to behave much like the KIM. In the phase diagram, a disordered 
high temperature phase and a low temperature phase with ergodicity breaking are separated by a critical line. The 
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phase transition is characterized by the same critical exponents 7] of the Ising model. After a quench into the low 
temperature phase, the non-stationary dynamics is analogous || to phase-ordering in the KIM. 

Despite these strong similarities, lack of detailed balance makes the IWDB in principle much different from the KIM 
and gives rise to new, interesting features that can be enlightened by the fluctuation dissipation relation. In the KIM, 
stationary states are equilibrium states characterized by time reversal invariance (TRI). Instead, in the stationary 
states of the IWDB, which are the counterparts of the KIM equilibrium states, TRI is violated and FDT (J2J) breaks 
down. The relation between xi^t^w) and C(t,t w ) may not be meaningful, as in systems with detailed balance, and 
—dx(C)/dC cannot be straightforwardly interpreted as a thermodynamic temperature. Interestingly, however, we 
can derive fluctuation-dissipation relations that generalize to the present model what is known in systems where 
detailed balance holds. In doing that, we uncover that the response function is not naturally related to the spin 
autocorrelation function C(t,t w ) — (&i(t)<Ji(t w )) , but rather to the quantity A(t,t w ) = {<Ji{t)(Ji(t w )(3i(t w )), Pi(t w ) 
being the space-time dependent inverse temperature. In stationary states the fluctuation-dissipation relation reads 

X (r) = A(0) - A(r) - A(t), (4) 

where A(r) is a quantity related to the lack of TRI, analogous to the asymmetry j^J in systems with detailed balance 
out of equilibrium. In the special case when the temperature is constant, (3i(t) = (3, the model reduces to the KIM. 
Since now TRI is recovered, one has A(r) =0 and then 

X (r) = A(0) - A(r). (5) 

Recalling that in this case A(t) — (3C(t) one recognizes the FDT, Eqs. (11121) . In the case with non-constant tem- 
perature, since A(t,t w ) ^ (3C(t,t w ) both in stationary and non-stationary states, the relation between x{i^w) and 
C(t,t w ) remains unclear. 

In the phase-ordering process following a temperature quench, two time quantities can be splitted into a stationary 
and an aging term. In particular, for the integrated response function one has x{^^w) = Xst(t — t w ) + Xag{t,t w ), 
where the aging part obeys the scaling form 

Xag (t,t w ) =t- a *f(y), (6) 

with y = t/t w , as generally expected [Tl| in phase-ordering systems. When detailed balance holds, the exponent a x 
is uniquely determined by the space dimensionality d and the dynamic exponent z. Since in the IWDB the value 
of z is the same Q of the KIM, one expects the same value of a x in the two models. In fact, in d = 2 we find 
results consistent with a x — 1/4, as in the KIM. This result complements those of Refs. 0,11] where it was found that 
the equilibrium critical exponents and the non- equilibrium persistence exponent of the IWDB were the same, within 
statistical errors, of those of the Ising model. The numeric results of this paper, therefore, strengthens the idea that 
the two models belong to the same non-equilibrium universality class. 

Despite this, the parametric plot of x(C) is different from that of the KIM. Interestingly, the shape of this function 
is similar to that found ^(J in a soluble model of sheared binary systems where detailed balance is also violated but 
local equilibrium still holds. In particular, the extrapolation of the numerical data to the large t w limit is consistent 
with the horizontal line typical of phase-ordering in the region of small C of the plot. 

This paper is organized as follows: In Sec. [H] we introduce the model. In Sec. IIIII wc derive a relation between the 
response function and correlation functions which generalizes the FDT to the present model. This relation allows to 
compute numerically the response function without applying a perturbation and to discuss the fluctuation-dissipation 
relation in Sec. IIVI In particular, the stationary states at high and low temperature are discussed in Sees. IIV A II 
IIV A 21 while the aging dynamics following a quench is considered in Sec. IIVBI A summary and the conclusions are 
drawn in Sec. IVl 

II. THE MODEL 

We consider the Ising model defined by the Hamiltonian 

H[a] = -J Wi = - H aiHi M ' ( 7 ) 

<ij> i 

where ai = ±1 is a spin variable on a d-dimensional lattice and < ij > denotes nearest neighbors i,j sites. Hi[a] = 
J^2j. crj i , where ji runs over the nearest neighbors of i, is the local field. 
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A dynamics is introduced by randomly choosing a single spin on site i and updating it in an elementary time step 
with a transition rate w([a] — ► [&']). Here [a] and [a 1 ] are the spin configurations before and after the move, which 
differ only by the value of <7j. In the IWDB transition rates are generic but their ratio must fulfill the condition 



w([a 



(8) 



With a constant 0%[a] = j3 one recovers the kinetic Ising model (KIM) in contact with a reservoir at the temperature 
T = /3 _1 . In this case Eq. (JSJ is the detailed balance condition with respect to the Hamiltonian J7J). In fact, one has 

- Y,P(°i H M - °'i H iW]) = P{H[a] - H[a']) (9) 

i 

and hence 

wQv] -> [a']) _ expi-pHP]} _ P«W\ (10) 



w([a']^[a}) exp{-0H[a]} P eg [a] ' 
where 



P eq [a](xexp{-l3H[a}} (11) 

is the canonical equilibrium probability. Detailed balance implies that stationary states of the model are also equilib- 
rium states with measure <|11[) . which are characterized by TRI. For a generic two-time quantity F(t, t w ), TRI implies 
F(t,t w ) = F(—t, —t w ). If TTI is also obeyed, namely F(t,t w ) = F(t — t w ), by shifting time by an amount t w , one 
also has F(t,t w ) — F(\ t — t w |). 

In this paper we consider the case when 0i[a] is not constant but depends on the configuration through the local 
field Hi[a], @i[o] = j3{Hi[a\). Physically, one can imagine a system in contact with reservoirs at different temperatures 
each of which couples to the spins with the same local field -Hj[er]. Notice that flipping ct; does not change Hi[a], 
hence (3i[<j) = f3i[<j')- However, this is true only for the site i where the flip occurs, while, in general, /3j[a] ^ 
for j ^ i. The very basic feature that makes this model different from the KIM is the fact that its transition rates 
do not obey detailed balance. This is expected on physical grounds, since different local temperatures in the system 
determine heat fluxes that break TRI and hence detailed balance. Mathematically, this happens because the argument 
of the exponential in Eq. © cannot be written as a difference H[a'] — H[cr], H being a generic function, as in Eq. JjJ), 
due to the factor /3i[<r]. Indeed, the term (3i[a]aiHi[a] is a function H[<r] of the configuration [a] but the term 
Si /3i[c]cr--ffi[°' / ] 7~t\ a '\ because it depends on both the configurations [a] and [a 1 ]. Because detailed balance is not 
obeyed, the stationary states of the model Q are not equilibrium states and, in principle, TRI is not expected. 

As discussed in 0,1a El, the present model contains, as special cases corresponding to particular choices of /3(\Hi\), 
the Voter, majority Voter and noisy Voter model, besides, clearly, the KIM. 



III. FLUCTUATION-DISSIPATION RELATIONS 



In this Section we derive a relation between the response function and particular correlation functions which 
generalize the result of Ref. 01 to the case of a non-constant (3i [a] . The derivation closely follows that of Ref . 0] to 
which we refer for further details. 

Let us consider a perturbing magnetic field switched on the j-th site in the time interval [t' , t' + At], 

hi(t) = h5 itj 6{t - t')9(t' +At-t) (12) 
where 8 is the Heavyside step function. The Hamiltonian J7J) is changed to 

H[a] =-jJ2 a W - Yl hi ^ ai = - £ °i H i[<r\> ( 13 ) 

<ij> i i 

where H^[a\ = JJ2j a U + hi(t). In the limit of vanishing h, the effect of the perturbation H12|) on the spin on site i 
at the time t > if is given by the linear response function 0, 

(14) 

h=0 



R lj{t ,t<)= Urn 

• n ' At^o At dhAt') 
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where here and in the following (. . .) means ensemble averages, namely taken over different initial conditions and 
thermal histories. Introducing the probability p([er],i) to find the system in the configuration [a] at time t, and 
the conditional probability p([cr],i|[er'],t') to find the configuration [a] at time t given that the system was in the 
configuration [a'] at t', the r.h.s. of Eq. (|I4|I can be written as 



d{a i {t)) 



(Tip([(7],t\[a'],t' + At) 



dp h ([a'],t' + At\[a"],t') 



dhj 



p([a"],t'). (15) 

h=0 



Here p and p h refer to the conditional probabilities of the unperturbed and perturbed system, respectively. Let us 
concentrate on the factor containing p h . The conditional probability for At sufficiently small is given by 

p h {[a'],t' + At\[a"],t') = + ™ h {[<*"] - M)Ai + 0(At 2 ), (16) 

where we have used the boundary condition p([cr'], t\ [&"], t) = S\ er n! (T /n. Furthermore, also the perturbed transition 
rates w h must verify the condition (JBJ), namely 

^pr^j] - «p {-£ AM^M - • ( 17 ) 

Expanding the perturbed transition rates in powers of h, one finds that the following form 

w h ([a'\ - [a"}) = w([a'} - [a"]) jl - - a'Jhj) } , (18) 

where iu([cr'] — » [a"]) are generic unperturbed transition rates obeying (JSJ, is compatible to first order in h with the 
condition (|17l) . 

Using Eqs. (|16[) and l|18|) . following Ref. the response function can be written as the sum of two contributions 

Ri,j{t,t') = [D itj {t,t',At) + A,,-(t,t',Af)] , (19) 

where 

A,i(M',At) = i £ - iP ([ -] > t|M,t' + At) £ <M!4«') (20) 



2 KM k"¥K] 



and 



Di,j(t,t',At) = \ J2 ^([ajM^t' + m^-^WjWMW''} - WMW'U). (21) 

[<r],K],[a"]#K] 

Using the time translational invariance (TTI) of the conditional probability p([cr], t\ [a'], t' + At) — p([a],t — At\ [cr'],t'), 
one can write Dij(t, t' , At) in the form of a correlation function 

D id (t,t',At) = -^(t-Ai)^')}, (22) 

where 

Bj =-E( ff i-^ikWH - KD- (23) 



Using Eq. (|TT)|l . Dij(t,t' , At) can be written as 



7),. ;( /.,'.A/,=i^j^ ,24, 



where 

AA itj {t,t') = + ^) - ^(t')]> (25) 
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Therefore, putting together Eqs. (|22|l . (|24() and taking the limit At — > we obtain 

RiAt,t') = l dAl ^ t,] - l«t)BAt% (26) 

where 

A 1J (M') = (ft(('K(%(0) (27) 
In the following, we will be interested in the integrated response function 

Xi, j ( t ^) = J t Ri,j(t,t')dt' (28) 

which correspond to the application of a perturbing field between the times t w and t. This quantity is easier to 
measure because switching on the perturbation for a finite time increases the signal to noise ratio. From Eq. I|26|) we 
have 

XiAt,t™) = \[Aij(t,t) -A itj (t,t w )} - i jf (ai(t)Bjtf)). (29) 

Eqs. 1|26I29[) are the principal results of this Section. They are relations between the response function and correlation 
functions of the unperturbed kinetics, which generalize the FDT. These relations hold both in stationary and non- 
stationary states, and do not depend on the choice of the unperturbed transition rates, provided the condition JHJ is 
obeyed. 

From Eq. (|29|l the integrated autoresponse function \(t,t w ) = Xm(^™)' which does not depend on i due to space 
translation invariance, is given by 



1 1 



X (t,t w ) = ~[A(t,t) - A(t,t w )] - - ^ (a^Bdt')), (30) 

where A(t, t w ) — Ai,i(t, t w ). We will use Eq. Q30|l for numerical computations in the next sections. As discussed in 0] 
this method to compute x(tyt w ) is much more efficient than traditional methods where the perturbation is switched 
on. 

In stationary states, a simplified expression for xi^^w) can be obtained which makes the role of TRI evident. In 
order to do this let us consider the integral 

1 r* 

I(t,tw)=2j dt (Bi(t)ai(t')) (31) 
Enforcing Eq. (|23|l . proceeding as in 01 the integrand can be written as 

( B, (iMt0) . "'"-y-') (32) 

Using Eq. (|32(l . replacing d/dt with —d/dt', due to TTI, and carrying out the integration one has 

I(t,t w ) = ~[A(t,t) - (AW^W^(^))] (33) 
Adding and subtracting I(t, t w ) on the r.h.s., Eq. (|30|l can be cast in the form |0J, with 

A(r) = 1 - (f3i(t w )ai(t)ai(t w )) + jf dt' [{a^B^t')) - (B^atf))] | . (34) 

If TRI is also obeyed, as in equilibrium states, one has (f3i(t)ai(t)<Ji(t w )) — {Pi{t w )<Ji(t)(Ji(t w )) and {ai(t)Bi(t')) = 
(Bi(t)(Ti(t')) , so that A(r) = 0. Eq. (0J becomes a linear relation formally identical to Eq. JSJ. When TRI does not 
hold, instead, A(r) ^ and the relation between x(r) and A(t) is no longer linear. As we will see in Sec. ljIV A|l . this 
is an efficient tool to check if a stationary state is invariant under time reversal and, if not, to quantify TRI violations. 
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FIG. 1: The phase diagram of the IWDB model where ti = tanh(2/3i) . The three broken curves correspond respectively, V to 
the voter model, I to KIM and M to the extreme or majority model. The two asterisks correspond to the two sets of parameters 
used in simulations. 



IV. NUMERICAL RESULTS 



In this section we present a numerical investigation of the dynamical properties of the model and, in particular, of 
the fluctuation dissipation relation (JSUJ). We chose unperturbed transition rates of the Metropolis type for single spin 
flip on site i 

«KM -> M) =mm{l,exp{-#M(ff[</] - H[a})}} , (35) 

which, as can be easily checked, obey Eq. Up down symmetry implies that Pi[a] does not depend on the sign of 
the Weiss field, (3i[cr] — (3(\Hi\). In the following we will consider a system on a square lattice in two dimensions. In 
this case the only possible values of the local field are Hi[a}/ J = 0, ±2, ±4; the model is then fully defined by assigning 
the three parameters /3(0),/3(2J),/3(4J). Moreover, with the Hamiltonian Q, the transition rates l|35(l do not depend 
on /3(0). Then, at this level, the couple of inverse temperatures (3(2 J), (3(4 J) is sufficient to characterize the model. 
However, j3(0) becomes relevant if the system is perturbed by an external magnetic field in order to measure response 
functions. Actually this quantity enters, through A(t,t') and Bi, the expressions (I26I29[I of the response functions. 
Then, response functions depend on (3(0), as already found numerically in Ref. [13j. 

The phase diagram of the IWDB was studied in Ref. 0,0- R was shown that in the plane of the parameters (3(2), (3(A) 
one can identify two regions separated by a critical line /3(4) = (3 C [(3(2)], as shown in Fig. QJ. The critical line starts at 
the inverse temperatures (3(2) = (l/2)atanh(l/2), (3(4) = oo, corresponding to the voter model, where the transition 
occurs in the absence of bulk noise, passes through the Onsager critical point with (3(2) = (3(A) = (l/2)arcsinh(l) 
and ends at (3(2) = oo, (3(A) ~ 0.22, corresponding to the extreme model, where the transition occurs in the absence 
of interfacial noise. For (3(A) < (3 C [(3(2)] one has an high temperature phase similar to the paramagnetic phase of the 
KIM. Starting from any initial state the system quickly attains a stationary state where the magnetization 

m = (o-j), (36) 

vanishes. For (3(A) > (3 C [(3(2)] there is a low temperature phase similar to a ferromagnetic phase. Here there are two 
possible dynamical situations, depending on whether the systems enters a state with broken symmetry, namely with 
m = ±Mbs 7^ 0, or not. In the former case a stationary state is entered; in the latter there is a phase-ordering process 
and the system ages. We will consider these cases separately in Sees. IIV AIIV Bl where we will present the results of 
numerical simulations of a two-dimensional system on a square lattice of size 1000 2 , with J = 1. 
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A. Stationary states 

1. 0(4) </?c|/9(2)] 

We have prepared the system in the stationary state at the inverse temperatures ,3(0) = 0.80, /3(2) = 0.44, /3(4) = 
0.30, which correspond to the paramagnetic phase. This state is quickly entered by the system by letting it to evolve 
from any initial condition. In the stationary state we checked that m = and that two-time quantities are functions 
of the time difference r alone. In the following, time will be measured in montecarlo steps. C(r), A(t) and x( r ) are 
shown in Fig. 

The behavior of C(r) is analogous to that observed in the KIM. Starting from C(0) = (af(t w )) — 1 the corre- 
lation function exponentially decays to zero. This is due to the decorrelation of the spin for large time differences 
lim r _ >00 C(t) = \\m T ^ 00 {ai{t w + T)ai(t w )) = (<r)(a) = m 2 = 0. Here we have introduced the simplified notation 
(a) = (&i{t)) to indicate that {<Ji(t)) does not depend on time t nor on site i due to TTI and space homogeneity. We 
will use this notation also in the following, dropping time and/or space variables whenever ensemble averages do not 
depend on them. 

A(t) behaves similarly. From the definition l|27|) . its equal time value is the average inverse temperature of the 
bath, A(0) = ((3) = 0.46. For large times difference also this correlation function decays to zero, since linv-xx, A(t) = 
m(a(3) = 0. Notice that C(r) and A(t) are proportional for large r but not for small r. This fact will be relevant 
in the following, when discussing the possibility to define a thermodynamic temperature from the parametric plots 
x(C), x(A) obtained plotting x{ T ) versus C(r) or A(t), respectively. 

The behavior of x( T ) is a l so closely related to what is known for the KIM. This quantity starts from x(0) = and 
saturates exponentially to a constant value Xoc- I n the KIM this value is the equilibrium susceptibility, namely the 
inverse temperature, Xqc = P- O ne could conjecture that in the IWBD this result can be generalized to Xoc = (/?)• 
However, as we will see shortly, this is not true due to lack of TRI. In order to discuss this point let us consider, in 
Fig. the parametric plot of x(A) . 

For the largest values of A, a relation formally identical to Eq. (J5J is obeyed. Recalling Eq. (0J this implies that the 
term A(r) is negligible. This, in turn, shows that TRI is satisfied in this time domain. This is reminiscent of what 
happens in out of equilibrium systems in contact with a single reservoir, where the linear FDT relation (j3J) is found 
on the r.h.s. of the x(A) parametric plot, despite the system is not in equilibrium. Since in this case A(0) = (f3) 
Eq. JSJ can be written as 

X (A) = (13) - A. (37) 

This shows that the average bath temperature (/?) enters the relation between x( T ) ancl A(t), as a natural general- 
ization of what happens in the KIM where one has x(A) = /3 — A. This behavior can be explained recalling that, in 
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A 

FIG. 3: The parametric plot x(A)- m the inset the small A sector is magnified. 



this sector of the plot, namely for r ~ 0, the response is provided by the fastest dynamical features. These are the 
microscopic flipping of single spins that locally and instantaneously equilibrate at the current inverse temperature 
(3i[cr\. Since the system is translational invariant, by taking ensemble averages one gets a sort of FDT with respect 
to the average bath temperature (/?), namely Eq. I|37[) . As larger values of r are considered, corresponding to lower 
values of A(r), slower dynamical features are probed which cannot follow the variations of /3i[<r]. In this sector, the 
contribution of A(t) becomes relevant and the parametric plot x(A) deviates from the straight line. Recalling the 
discussion of Sec. A(t) ^ is related to the lack of TRI. These considerations, then, suggest that the paramet- 

ric plot can be used as a convenient tool to detect and quantify TRI violations in non equilibrium stationary states. 
Breakdown of TRI is also responsible for the saturation of x(r) to a value Xoo < (/?)■ 

The calculations of Sec. IIHI clearly show that A(t,t w ) is the correlation function naturally associated to x(t,t w ), 
rather then the autocorrelation function C(t, t w ) which does not enter the generalization of the fluctuation dissipation 
theorem <|30|) . Nevertheless, we consider, in Fig. also the parametric plot of x(r) versus C(r), since, by analogy 
with systems with detailed balance, this relation is often considered in the literature |13| . 

Also in this case, for the larger values of C, the curve x(C) obeys a linear relation x(C) = (3 ■ (1 — C), with (3 = 0.55. 
Although this fact has suggested the interpretation [RJ of as a thermodynamic temperature, since C(t,t w ) does 
not enter the fluctuation dissipation relation (|30[) . this reading is unmotivated. Notice in fact that (3 ^ ((3). Indeed, 
one should have (3 — (f3) if C(t) oc A(t) in the regime considered, namely for small r. Instead, as discussed previously, 
this is not the case. Actually, we have checked that (3 ^ {(3) even when the choice of (3(0) proposed in 0] is made, 
both in the high and low temperature phases. 

For small values of C the curve x(C) strongly deviates from the straight line. The parametric plot x(C) can be 
compared to that found in stationary states of other systems without detailed balance and, in particular, in an exactly 
soluble model of binary systems under shear flow [lC| , where a similar pattern was found. 

2. /3(4) > /3c[/3(2)] 

We have prepared the system in the stationary state at the inverse temperatures (3(0) — 0.80, (3(2) — 0.68, (3(4) = 
0.37, corresponding to the ferromagnetic phase. This state is quickly entered by the system by letting it to evolve 
from any initial condition with a broken symmetry m ^ 0. We consider the stationary state where the magnetization 
attains the positive value m = Mbs — 0.83. Two time quantities, denoted by Cbs(t), Abs(t), Xbs(t) are plotted in 
Fig- El 

Cbs(t) behaves similarly to the KIM. It decays from Cbs(0) = 1 to the large r value Cbs(°°) — ( fJ )( cr ) = M BS — 
0.69. Analogously, A B s(t) decays from A BS (0) = ((3) = 0.45, to A BS (°°) = (°-)(cr/3) = M BS (ofl) = 0.30. xbs(t) 
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FIG. 4: The parametric plot x(C). 
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FIG. 5: C(t),A(t) and x( T ) are plotted against r. 



grows from zero up to the constant value Xoo • If TRI were obeyed, from Eq. @ one should have \oo = (/?) — Mbs (c/3) . 
Indeed, for the KIM this equation gives Xoc = — Mbs) which is m fact the equilibrium susceptibility. However, 
TRI is violated in the IWDB, as it is clear from the fact that the parametric plot of xbs(t) versus Abs(t), shown in 
Fig. El is not a straight line. 

Deviations from the straight line are due to the term A(t) in Eq. Q), which signals the breakdown of TRI, and 
makes xss(t) saturate to a value Xoo < — Mbs(&(3)- Notice however that, also in this case, for the largest values 
of Abs, the linear relation xbs(Abs) = (P) — Ass is obeyed, as in the paramagnetic phase, implying that TRI is 
satisfied in this time domain and that the average temperature can be extracted from this region of the plot, in 
analogy with the KIM. 

In Fig. [7] the parametric plot of Xbs( t ) versus Cbs( t ) is also shown. 
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FIG. 6: The parametric plot Xbs(Abs)- In the inset the small Abs sector is magnified. 




FIG. 7: The parametric plot Xbs(Cbs)- 



This plot is analogous to the one found in the paramagnetic phase, and similar considerations can be made. In 
particular, we hnd xbs{Cbs) = $ • (1 — Cbs): with (3 = 0.54, for the largest values of C. As already discussed in 
Sec. HIV A 10 . there is no reason to interpret this quantity as an inverse temperature, and again (3 ^ (/?). 

Let us also introduce the connected two time quantities, that will be used in Sec. lIVBl Using the general definition 
of the connected correlation function V between two observables O and O', V = {00') — (0)(0'), the connected 
two-time quantities associated to Cbs(t) and Abs(t) are 

C BS (r) = C BS (r) - M% s , (38) 



Abs(t) = A bs {t) - M bs (<t(3) 



(39) 
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B. Aging dynamics 

In this Section we study the non equilibrium process following a quench from an initial disordered state with m = 
to the final inverse temperatures /3(0) = 0.80, (3(2) = 0.68, /3(4) = 0.37. Notice that these are the same temperatures of 
the previous section HV A 21 corresponding to a point in the orderedphase. In this case one observes a phase ordering 
process where domains of two phases with m = ±Mbs coarsen In the interior of such domains the system is 
found in the stationary state studied in the previous section. In analogy to what in known for the KIM, and more 
generally in aging systems |Tll |. we expect quantities such as the equal time correlation function 

G{r,t) = {a i {t)a j {t)), (40) 

i and j being two sites whose distance is r, or two-time correlation functions and response to take the additive 
structure [12| 

G(r,t)=G st (r) + G ag (r,t), (41) 

C(t,t w ) = C st (T) + C ag (t,t w ), (42) 

A{t,t w )=A st (r)+A ag (t,t w ), (43) 

X(t,t w )=Xst(T)+Xag(t,t w ). (44) 

The presence of the stationary state in the bulk of the growing domains is the origin of the contributions 
G st (r),C st (T), A st (T),Xst(r), while the terms G ag (r,t),C ag (t,t w ), A ag (t,t w ),x ag (t,t w ) take into account the aging 
degrees of freedom in the system. 

In the KIM quenched to the final temperature G s t(r) is the correlation function of the stationary state with 
broken symmetry at the same temperature /3 , namely the correlation Gss(r). We define G s t(r) in complete analogy 
for the IWDB, GbsI 7 ") being the quantity measured in the stationary state at the same inverse temperatures. G ag (r, t) 
can then be obtained by subtraction, by using Eq. (|41l) . In the scaling regime G ag (r,t) obeys 

G ag (r,t) = M BS g(x), (45) 

x = XTtJ' This property will be tested below. 

The typical size of domains can then be computed as the half height width of G ag {r,t). This quantity is shown in 
Fig- [HI After the initial transient, L(t) has a power law behavior L(t) ~ t^.We measure 1/z = 0.50, as for the KIM. 

Coming back to the scaling l|45|) . in order to check this form we plot in Fig. [5] G ag (r, t)/M BS against x for different 
values of t. According to Eq. I|45|) one should find data collapse for different times. Actually, the collapse is good even 
if worst than in the KIM, particularly for x ~ 0. Notice also that the form of the scaling function g(x) is very similar 
to that of the KIM. 

Let us discuss now two-time quantities. Analogously to what discussed above for G(r, t), in the KIM one has 

C st (r) =C bs {t). (46) 

We define C s t(r) in complete analogy for the IWDB, and the same is assumed for Ast(r), 

A st {r) = A bs {t). (47) 

For the integrated autoresponse function, in systems with a constant /3, Xst(i~) is the response produced in the bulk 
of domains and is defined as the quantity related by Eq. |0J to the stationary parts of the correlation functions. In 
this case A(t) = 0, Eq. gj is the FDT © and 

Xstir) = A st (0) - A st (T) = A BS (0) - Abs(t) - A BS {0) - A bs {t). (48) 

Clearly, since x s *(t) is related by FDT JSJ to the correlation function of the broken symmetry equilibrium state, it is 
the integrated autoresponse function of that state. Then one has 



Xst(r) = Xbs(t)- 



(49) 
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FIG. 9: G ag (r,t) is plotted against x for different times t generated from t n =Int[exp(n/2)] with n ranging from 14 to 20. In 
the inset G ag (r,t) for the same values of t is plotted for the KIM quenched at T = 0. One observes the same small x behavior 
g(x) ~ 1/x. 



For the present model, in full analogy to the case with constant 0, we use Eq. 149|l to define Xst(r). Namely, Xst{i~) is 
the quantity measured in the stationary state at the same inverse temperatures in the previous section. Let us now 
turn to discuss the properties of the aging contributions in Eqs. H42I43I44[) . Recalling the behavior of Cbs{ t ) onc 
concludes that C ag (t,t w ) decays from C ag (t w ,t w ) = M% s to zero, as in the KIM. In analogy to the KIM we expect 
it to obey the scaling form 

C ag (t,t w ) = h c (y), (50) 

where y — t/t w , with the power law hc(y) ~ y~ x / z for large y, A being the Fisher-Huse exponent. Analogously, given 
the behavior of Abs( t ) discussed in the previous Section, one concludes that A ag (t, t w ) decays from A ag (t w , t w ) — 
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FIG. 10: C ag (t,t m ) is plotted against y. 



Mbs^P) = 0-30 to zero. We expect a scaling form 

A ag (t,t w ) = h A (y), (51) 

as for C ag (t,t w ). The response Xa g (t,t w ) is produced by the interface degrees of freedom whose number goes to zero 
during the ordering process. For this reason, in the d = 2 KIM this quantity after reaching a maximum for y ~ 1 
decays to zero. A scaling behavior is obeyed, namely 

Xag(t,t w )=t- a xf(y), (52) 

with the power law f(y) ~ y~°" x for large y and a x consistent with a x = 1/4 [l8j. In order to check these scaling 
forms we have extracted the aging terms as 

C ag (t,t w ) = C(t,t w )-C st (T)=C(t,t w )-C BS (T), (53) 



A ag (t,t w ) = A(t,t w ) - Astir) = A(t,t w ) - A BS (t), (54) 

and 

Xa 9 (t,t w ) = x(t,t w ) - Xst(r) = x(t,t w ) - Xbs(t). (55) 

According to Eq. 1)50(1. curves C ag (t, t w ) corresponding to different values of t w should collapse when plotted against 
y. This type of plot is shown in Fig. QJJI 

This figure shows that the data collapse is not very good for all times . This can be associated to the presence of 
preasymptotic effects. However the quality of the collapse gets better for the largest values of t w and t. Indeed, while 
the two curves with the smallest values of t w (t w = 100, 200) do not collapse at all, there is a tendency to a better 
collapse as t w gets larger. For the two largest values (t w = 1600, 3200) one has a nice collapse from y ~ 5 onwards. 

A similar situation is observed for the correlation A ag (t,t w ), as shown in Fig. 1111 Similarly to what found in 
stationary states, we find that C ag (t,t w ) cx A ag (t,t w ) for large y. 

According to Eq. (|52*|) . the exponent a x can be extracted as the slope of the double-logarithmic plot of Xag(t,t w ) 
against t w , by keeping y fixed. We do this in Fig. 1121 for different choices of y. We observe a good power law behavior, 
for every value of y. Best fit exponents are in the range [0.23 — 0.28], depending on y, suggesting that the same value 
a x — 1/4 of the KIM is found here. Then, in order to check the data collapse, in Fig.^|we plot t^Xag(t,t w ), with 
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w 

FIG. 12: Xag{t, tw) is plotted against t w for fixed y — t/t w , and y — 3,5, 7, 9, 11, 13, 15, 17, 19 from top to bottom. The dashed 
_ i 

line is the power law t w 4 . 



a x = 1/4, against y, for different t w . The collapse is indeed rather good for the two largest values of t w , implying 
that Eq. I|52l) with an exponent consistent with a x — 1/4, as in the KIM, is asymptotically obeyed. This result 
complements those of Refs.0,@ where it was shown that the IWDB has the same equilibrium critical exponent and 
the same persistence exponent of the KIM. This strongly indicates that this two model belong to the same equilibrium 
and non- equilibrium universality class. 

In Fig. El the parametric plot of x(^) is shown. 

In order to understand this plot one has to consider separately the short time separation regime (ST), namely the 
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FIG. 13: twXag{t,t w ) is plotted against y — t/tw The dashed line is the power law y 




limit t w — > oo with r/t w <C 1, and the large time separation regime (LT), where t w — * oo with t /t w ^> 1. In the 
ST, given the scalings I|5UI51I52I) the aging parts of two time functions remain equal to their equal time value. For 
A ag (t,t w ) one has A ag (t,t w ) = A ag (t w ,t w ) = M B s{^P)- Then, from Eq. g3J| one has A(t,t w ) = M BS ((jf3) + A^(t) = 
Mbs{&P) + Abs{ t )- Recalling the behavior of Abs{^i tw) one concludes that, in the ST, A(t,t w ) decays from 
(j3) = 0.45 to Mbs{&(3) = 0.3. In this time domain one has Xag(t,t w ) — and hence x(t, i w ) = Xs*( r ) = Xbs(t). 
Therefore, on the r.h.s. of the parametric plot of Fig. El for A > 0.3, one should find exactly the same curve found 
in the stationary broken symmetry state, namely Fig. This curve is the broken line in Fig. 1141 This implies that 
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FIG. 15: The parametric plot x(C)- The dashed line is the expected asymptotic behavior, for t w = oo. 



for the largest values of A, let's say for A > 0.4, Eqs. I|5I37|I are obeyed, as discussed in Sec. IIV A~2l 

The numerical simulations can only access finite t w , and some deviations from the asymptotic curve are then 
observed. Notice, in particular, that the equal time value A(t w , t w ) has a weak time dependence. In this case there is 
a monotonous decrease of this quantity. Recalling that A(t w ,t w ) — ((5(t w )) this means that in the kinetic process the 
average temperature is slightly increasing. This happens because in the phase-ordering process the fraction of bulk 
spins, with /Jj = (5(4), grows in time. Since in this case /3(4) < (5(2) < (5(0) this corresponds to an increase of the 
average temperature. Despite these finite time effects however, the data clearly show that the curves converge to the 
broken line increasing t w . 

On the left hand side of the plot the LT is probed. In this regime Ass(t, tw) — 0, and A(t, t w ) = A ag (t, t w ) decays 
from Mbs(&I3) to zero. On the other hand Xst(r) has already reached its asymptotic value Xst( T ) — Xoo so that 
x(t,t w ) = Xoo + Xag(t,t w ). According to the scaling form (|52|l . Xag(t,t w ) vanishes in the large t w limit. Then, for 
t w = oo on the left hand side of Fig. [21 one should hnd the horizontal straight line typical of phase-ordering systems. 
For finite values of t w , Xag(t,t w ) still contributes to the response and the curve overshoots the asymptotic value Xoo- 
However, as shown in Fig. 1141 the asymptotic curve is approached increasing t w . 

Let us now consider the parametric plot of x(£, t w ) versus C(t,t w ), shown in Fig. 1151 

Repeating the same considerations as for the previous figure, one concludes that on the r.h.s. of the figure the 
curves approach the curve of xbs( t ) against Cbs( t ) of Sec. II V A"21 namely Fig. [7] in the large t w limit. We stress 
again that, as already discussed in Sec. IIV Al there is no reason to associate the quantity extracted from this 
sector of the plot, to a thermodynamic temperature, as claimed in |l3j . 

On the l.h.s. of the figure our data are consistent with a convergence to the flat line %(C) = Xoo, typical of 
phase-ordering. Notice, however, that the whole shape of the parametric plot is different from that of the KIM, due 
to the different relation between the stationary parts of x(~t,t w ) and C(t,t w ), which shows up in the large-C region. 
Moreover, also in this case the plot is similar to what observed in binary systems under shear [lCJ. 

V. SUMMARY AND CONCLUSIONS 

In this paper we have studied a modified Ising model, the IWDB, where the temperature entering the transition rates 
depends on space and time through the system configuration and detailed balance is violated. This model is known 
to share many properties of the Ising model, including the phase diagram, critical exponents 7] and non-stationary 
dynamics |8|]. 

In systems with detailed balance a relation between the integrated response function x(t,tw), the autocorrelation 
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function C(t,t w ) and the asymmetry A(t,t w ), a term related to the possible lack of TRI, can be obtained under 
general assumptions jot Il4|. This fluctuation-dissipation relation applies also in non-stationary states, namely out of 
equilibrium. 

In this paper we have derived an analogous fluctuation-dissipation relation for the IWDB. The result is similar to 
the case with detailed balance but the role played by C(t,t w ) is now played by the correlation A(t,t w ) between the 
spins and the time dependent local inverse temperature. Since (3i{t w ) enters the transition rates and is, therefore, 
correlated to the spin configuration, A(t,t w ) is not simply related to C(t,t w ). It is therefore natural to consider the 
relation xiA) f° r which a generalization of what is known in systems with detailed balance seems to be possible, 
instead of the relation x(C) whose meaning, as far as we can see, remains unclear. 

In the stationary states of the model, which are the counterparts of the Ising equilibrium states, the fluctuation- 
dissipation relation (0} is formally similar to the FDT JSJ, with the important difference of a non vanishing A(t, t w ), 
determined by the violation of TRI. This term makes x{A) non linear. However, for small time differences, namely for 
the largest values of A, the asymmetry can be neglected and one recovers a linear relation, as in equilibrium systems, 
with the average inverse temperature {(3) playing the role of f3 in equilibrium systems. 

After quenching the systems into the ferromagnetic phase, a non stationary process is observed, similar to the phase- 
ordering kinetics of the KIM. We find that the response function exponent takes a value consistent with a x = 1/4, 
as for the KIM. This fact complements previous results @, @ on the universality between the two models, both in 
equilibrium and out of equilibrium. The shape of the plots x(^) an( i xiP) can a l so be discussed in strict analogy 
to what observed in the KIM. In particular, one finds the flat horizontal line typical of phase ordering systems. 
Interestingly, the parametric plot x(C) is similar to that of a soluble model of sheared binary systems where detailed 
balance is also violated, suggesting some qualitative similarity between these two model. 
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